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Abstract 

Any quantum computational network can be constructed with a sequence 
of the two-qubit diagonal quantum gates and one-qubit gates in two-state 
quantum systems. The universal construction of these quantum gates in the 
quantum systems and of the quantum computational networks with these 
gates may be achieved with the help of the operator algebra structure of 
Hamiltonians of the systems and the properties of the multiple-quantum 
operator algebra subspaces of the Liouville operator space and the specific 
properties of the quantum algorithm corresponding to the quantum network. 
As an example, the two-qubit diagonal gates are exactly prepared in detail 
in superconducting Josephson junctions. 



Quantum computation becomes increasingly attractive largely due to the 
reason that it can be much more powerful than the classical counterpart and 
has an extensively potential application. 1-4 Quantum computational net- 
works are usually composed of a sequence of quantum gates that could be 
built up with one- and two-body interactions of quantum systems. 5-7 The 
practical quantum gates should satisfy several requirements. For example, 
they should be simple so that any quantum computational network can be 
built easily out of these gates 5> 8 and could be conveniently prepared experi- 
mentally in those accessible quantum systems. Owing to inevitable decoher- 
ence and dephase effects in practical quantum systems the quantum gates 
should be chosen suitably in order to minimize thses effects on the quantum 
computation built up with these gates. The practical quantum gates may in- 
clude the three-qubit Toffoli gate, 9 the Deutsch there-qubit universal gate, 5 
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the two-qubit universal gate, 10,11 the two-qubit XOR gate 12 and the two- 
qubit diagonal gate, 13 etc. In recent years, a variety of two-state quantum sys- 
tems have been investigated experimentally to construct the quantum gates. 
These systems include trapped ions, 14 cavity quantum electrondynamics, 15 
nuclear spin systems in molecules, 16 ' 17 and the superconducting Josephson- 
junction arrays. 18-21 Quantum computation obeys quantum physical laws 
and time evolution of a quantum system during quantum computing is de- 
scribed by the Schrodinger equation, 22 while a quantum algorithm whose 
quantum network usually consists of a sequence of unitary transformations 
places constraints on the time evolution and hence constraints the form of the 
effective Hamiltonian of the quantum system in the quantum computing. 13 A 
quantum algorithm is usually more dependent on mathematical and quantum 
physical principles and largely independent of a general quantum system, but 
the quantum-gate construction is always performed in a practical quantum 
system. Therefore, in real quantum computation it is important to consider 
the specific properties of both quantum systems and quantum algorithms for 
the unified and systematical construction of the quantum gates in a variety of 
accessible quantum systems and of any quantum networks with these gates. 
However, this point is often neglected and is really not paid much attention 
to. In fact, such a universal construction may be achieved simply with the 
help of the operator algebra structure of the Hamiltonian and the properties 
of the multiple-quantum operator algebra subspaces of the Liouville operator 
space of a quantum system, and the specific properties of a given quantum 
algorithm. 23 In this communication it is shown how to build exactly the two- 
qubit diagonal quantum gates in superconducting Josephson junction arrays 
18-21 ma ki n g use Q f the operator algebra structure of Hamiltonian and 
the properties of the multiple-quantum operator algebra subspaces of the 
quantum system. This general principle allows one to apply to other coupled 
two-state multiparticle quantum systems. 

The quantum computational network of a given quantum algorithm usu- 
ally consists of a sequence of quantum circuit units. 5 Each such unit (e.g., 
the kth unit) is described by a unitary operator that may be written as an 
exponential form: = exp(-iHktk), where is the time-independent 

effective Hamiltonian associated with the fcth circuit unit. To implement the 
quantum circuit units in a quantum system one possible strategy is that the 
unitary exponential operator Uk(tk) is first decomposed completely into an 
ordered product of a series of elementary building blocks G m (X m ), i.e., the 
one-qubit gates and the two-qubit diagonal quantum gates 13 



2 



U(t) =n G m (\ m ) (A) 

m 

The two-qubit diagonal gate G k i(\ki) is denned in an operator form by 

Gki{Ki) = exp(-i\ kl 2I kz I lz ) (Bl) 
or in the unitary representation in the conventional computational basis by 
G k i(X k i) = Diag[e~^ Xkl , e^ A «, e ^ Xkl , e~^ Xkl \ (B2) 
where the operators I pz = \o vz (p — k, I), a is Pauli s operator and the two- 
qubit diagonal gate acts on only the fcth and Zth qubits in the quantum 
system simultaneously. The quantum computation then can be carried out 
by implemeting the decomposed unitary operations U k (t k ) of Eq.(A) in a 
feasible quantum system. The decomposition of Eq. (A) may be achieved in 
an exact and unified form and may be simplified greatly with the help of the 
operator algebra structure of the Hamiltonian Hk and the properties of the 
Liouville operator space and its multiple-quantum operator subspaces. 23 If 
the Hamiltonian Hk is a member of the /ongitudinal magnetization and spin 
order {LOMSO) operator subspace the decomposition of Eq.(A) is simple 
due to the fact that each pair of the base operators of the operator subspace 
commute. When the Hamiltonian H k is a zero-quantum operator, the de- 
composition is carried out in the zero-quantum operator subspace, while H k 
is an euen-order multiple- guantum (EVOMQ) operator the decomposition 
may be achieved in the EVOMQ operator subspace. For example, all the 
diagonal phase quantum gates can be constructed in the LOMSO operator 
subspace, 3-5,13 while an arbitrary SWAP operation 5 between any pair of 
quantum bits may be built up in the zero-quantum operator subspace. 

The preparation for the elementary building blocks, that is, the two-qubit 
diagonal gates in a variety of two-state coupled multiparticle quantum sys- 
tems may also be simplified with the help of the operator algebra structure of 
the Hamiltonians of the quantum systems, and the properties of the Liouville 
operator space and its multiple-quantum operator subspaces. For example, 
for a weak coupled N-spin (1=1/2) system 24,25 its Hamiltonian is a memeber 
of the LOMSO operator subspace. The two-qubit diagonal gate then can be 
achieved with the aid of the properties of the LOMSO subspace that any pair 
of base operators of the subspace commute, and may be prepared exactly by 
using a sequence of spin echo units with selective radiofrequency pulses in the 
system. 13 Likewise, if the Hamiltonian is an element of the zero-quantum op- 
erator subspace the elementary building blocks may be prepared by utilizing 
the properties of the zero-quantum operator subspace. The Hamiltonians of 
the superconducting Josephson junction system may consist of the multiple- 
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quantum operators (see below). Then the properties of the multiple-quantum 
operator subspaces, e.g., the EVOMQ subspace will be helpful for the exact 
preparation of the two-qubit diagonal gate in the system. 

It has been suggested in recent years that superconducting Cooper-pair 
boxes with Josephson junctions could be used as qubits. 18-21 Very recently, 
the coherent time evolution of the two charge states in a single-Cooper- 
pair box is observed experimentally. 21 This demonstrates for the first time a 
practical solid-state qubit for quantum computation. It is discussed in de- 
tail below how to prepare the two-qubit diagonal gate in the low-capacitance 
SQUID-controlled Josephson junction arrays. 18,20 Obviously, any one-qubit 
gate is easily implemented in the system and will not be further discussed. 
Now consider a system of two coupled SQUID-controlled Josephson junc- 
tions, without loss of generality. The Hamiltonian for the system can be 
written generally in a spin (1=1/2) language 18 ' 20 

H T = E k h (V Xk )a k + E[ h {V Xi y z - \E k ^ Xk )a k x - \E l j(<$> Xl )a x 

-[E k ^ Xk )E l j^ Xl )/E L )a k a y (1) 
where a is Pauli s operator. By a simple unitary transformation: H = 
UyHxUy where U y = exp(— i(f> k I ky ) exp(— ifaliy), the Hamiltonian (1) can 
reduce to a simpler form 

H = Vlkhz + tyhz + K.Jkl2IkyIly (2) 

where the operators l iiX = \a\ , tan& = 2E t ch ( y V Xt )/E i J (<S> Xi ), 
Qi = -2E l ch (V Xi ) sin0j - E l j($ Xi ) cosfa (i = k,l; fx = x,y,z), and nJ kl = 
—2E k (V Xk )Ej(V Xl )/E L . By employing the propagator U(t) = exp(-iHt) 
corresponding to the Hamiltonian (2) one can prepared the two-qubit diag- 
onal gates based on the Baker-Campbell-Hausdoff formula, correct to the 
third order approximation, 13 

Gki(Xki) = exp(-i|F x ) exp(-iHt/2) exp(-i7fE z ) exp(-iHt) 
x exp(i7rF^) exp(—iHt/2) exp(i^F x ) 
= eM-i27cJ kl 2hzIizt)+0(t 3 ) (3) 
where the operators F^ = I kfM + (/i = x, y, z) and the unitary operations 
exp(±iaF Al ) are single-qubit operations. 

Actually, the two-qubit diagonal gates of Eq.(B) can be prepared ex- 
actly with the help of the operator algebra structure of the Hamiltonian 
of the superconducting Josephson junction arrays and the properties of the 
multiple-quantum operator subspaces. By a simple transformation: H e = 
exp(i^F y )H exp(— i\F y ), the Hamiltonian (2) can be converted into a mem- 
ber of the EVOMQ operator subspace 
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H e — Qklkz + &lhz + ^Jkl{Ikxhx + Ikyhy) + KJ k i{ — 

Ikx hx Ikyhy) (4) 

where I kz and I\ z are the longitudinal magnetization operators, (Ikxhx + 
Ikyhy) and {—Ikxhx + Ikyhy) the zero- and double-quantum operators, respec- 
tively, 24,25 indicating that the operator H e is an even-order multiple-quantum 
operator. Then by utilizing the properties of the EVOMQ operator subspace 
one can diagonalize the operator H e in the operator subspace. 23 The result 
turns out to be as follows 

h = v + H e v = n' k i kz +n' l i lz (5) 

where the unitary operator V is given by 

V = exp(—iaQ ) exp(-if3Q 2 ) (6) 
with the zero-quantum operator Q and the double-quantum operator Q 2 : 

Qo — ^{Ikxhy — Ikyhx) and Q 2 = 2{I kx h y + Ikyhx)- 

For simplification, the ratio 7 between the parameters Q k and Q l in Eq. (5) 
is set as 7 = (2m + l)/(2m) (m = ±1, ±2, ...), then 

n' k = 7^;, (r) 

and in the system of the two coupled SQUID-controlled Josephson junctions 
the evolutional time t is chosen as t p so that 

tt' k t p = (2m + l)n (8a) 

and hence 

Sl'jp = 2mn. (86) 
Note that Q k and Q t are dependent on the parameters Q k , Qi, and J k i. If 
Qi (i = k, I) and J k i are adjusted suitably, this could be achieved by tuning 
the external flux $ Xi and the inductance L in the superconducting Josephson 
junction circuit, 18 ' 20 equation (7) can be met for a given ratio 7. By Eqs.(5)- 
(8) the propagator corresponding to the Hamiltonian H e of Eq.(4) can take 
a simpler form 

exp(-iH e t p ) = V 2 exp(-i-Kl kz ). (9) 
In order to prepare the two-qubit diagonal gate an echo sequence is con- 
structed as follows 

Pi = exp(-iH e t p ) exp(-iirl kx ) exp(-iH e t p ). (10) 
With the help of Eqs. (6) and (9) this operation is simplified 

Pi = exp[-i8(-a+/3) Ikyhx] exp(-inl kx ). (11) 
Then one obtains the two-qubit diagonal gate from Eqs. (10) and (11) 

G k i{\ k i) =Tiexp(-iHt p )exp(iirI kz )exp(-iHt p )Tf exp(iirl kx ) _ (12) 
where 7 = (2m + l)/(2m), the evolutional time t p = (2m + l)n/Q' k , the 
parameter X ki = 4(— a + /3), and the unitary operation Ti represents the 
single-qubit composite rotation operation: 

Ti =exp(-qi kz )exp(qi kx ). (13) 
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To determine the parameters Q' k , Q) and a, (3 one needs to expand the trans- 
formation V + H e V of Eq.(5) by using the rotation transformation between 
any two product operators. 25 The equations to determine these parameters 
are therefore obtained explicitly 

Qk sin ai cos a 2 + &i cos a\ sin a 2 — ^Jm sin on sin a 2 = (14a) 
flk cos ai sin a 2 + sin a.\ cos a 2 + tt^m cos a-y cos a 2 = (146) 
and 

Q k = flk cos cki cos a 2 — fi; sin cti sin a 2 — 7r cos a\ sin a 2 (15a) 
f2 z = Qi cos cti cos a 2 — sin cti sin a 2 — 7r J^^ sin a\ cos a 2 (156) 
where oi\ = —a + /3 and a 2 = a + (3. It is easy to find two independent 
solutions a 2 = af to Eqs.(14) that are explicitly expressed as 

tanaf = \(-8±^5 2 + A) (16a) 
and the corresponding a± — af that are determined by 

tanaf = - I^(-o± V^TI) (166) 

where 5 = [il 2 —ill+(TcJ k i) 2 ]/(7rJki^lk)- Note that the propagator exp(—iH e t p ) 
is the same for a given time t p for the two solutions (af , af) and (cuf, a 2 ). 
Without lossing generality, here the solution is taken as (af, af) and the 
parameters fl k and Q l are calculated by inserting the solution into Eqs.(15). 
Consequently, one can explicitly obtain equation (7) that constrains the pos- 
sible values of the physical parameters Q k , ^i, and J k i for a given 7 = 7+ (it 
can prove that 7 7^ 0, ±1). In theoretical design the parameters \ k i = Aaf 
and 7 + = (2m + l)/(2m) are given in advance. Then equations (7) and 
(16b) can be simplified and recast a one- variable (//) quadratic (the param- 
eters /i = fli/Qk and v = nJki/flk below): 

[p - tan(±A fc ,)] V - 7 + {l + 2p[p - tan(iA«)] + tan 2 (±A fci )}/i 

+ (j9 7 +) 2 + 1 = (17a) 

with 

= tan(\\ kl )[l +tan 2 (|A fc Q] 

P 2[l+tan(|A fc ,)tan(iA fc ,)] 
After the parameter /x is determined from Eq.(17a) another parameter v can 

be obtained by 

v = p{ri + -n) (176) 
In Eqs.(17) it is required that 

/i 7 + [l + tan 2 (±A H )] > 1 
and the solution (af, af) also requires that 
a 2 + is 2 - 1 

^ > 2//tan(±A w ) + 2z/. 
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Given the parameters Am and 7 + , the parameters fi and v can be determined 
by Eqs.(17) and consequently, the physical parameters Qk, Qi, and J^i in the 
superconducting Josephson junctions are set up properly, while the parame- 
ter fl k is calculated through Eq.(15a) by using these determined parameters. 
The evolutional time t p = (2m + l)n/Q' k is therefore obtained exactly. 

In practice, for a given set of the physical parameters {flk, -hi} ob- 
tained according to the above theoretical design, the Hamiltonians H and H e 
are determined by Eqs.(2) and (4), respectively. By choosing the evolutional 
time as t p = (2m + l)7r/f2' fc the propagator exp(—iH e t p ) will take really ei- 
ther the form: V 2 exp(— iirlkz) (Eq. (9)) or the form: V 2 exp(iirlk z )- Although 
the latter form is different from Eq.(9), the same result of the echo sequence 
Pi of Eq.(ll) is obtained. On the other hand, the parameter \ M in Eq.(12) 

will take really rnr + Aaf, where af = arctan{ — — — — (—5 + \J 8 2 + 4)} 

ill 2 i li 

is determined from Eq.(16b) by the above theoretical design, but the factor 
nir does not give a net effect on the two-qubit diagonal gate of Eq.(12). Thus, 
the two-qubit diagonal gate Gki(Aki) (^ki — ^ a t) is really prepared exactly 
by the echo sequence (12). 

As an example, to prepare the two-qubit XOR gate 12 whose unitary 
operation can be expressed as 

Uxor = exp(if) exp(-zf^) exp(-i|/ fe2 ) exp(-if J te ) 
x ex.p(-inl kz l lz ) exp(i|/ iy ) 
one needs to prepare the two-qubit diagonal gate with = |. For \m = f 
and 7 + = I (m = 1) the parameters ji and v are obtained from Eqs.(17): \i = 



and v = and Q' k from Eq.(15a): §^ = ±\ ^26 - by/2, respec- 

tively. The evolutional time is therefore determined by |f2J t v = , 127r 
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